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LOCALITY  PRINCIPLE  IN  WAVE-llECHANICS 
(Wave  Mechanics/Local  Density  of  States/Local  Partition  Function/ 
Cluster  Methods/Moment  Methods/Path  Integral  Methods) 

W.  Kohn  and  A.  Yaniv 

Department  of  Physics,  University  of  California  San  Diego, 

La  Jolla,  California  92093 


ABSTRACT 

This  paper  proves  a locality  principle  for  a wave-mechanical 
particle  governed  by  the  Schroedinger  equation.  It  is  shown 
that  p(r,g),  the  Laplace  transform  of  the  local  density  of 
states  n(r,E),  depends  significantly  only  on  the  potential 
V(r')  at  points  r'  near  r.  The  effect  of  changes  of  V(r') 
at  distant  points  r'  (|r'  - r | > d)  on  p(r,8),  decay  in  a 
Gaussian  fashion  with  a.  This  result  sheds  some  light  on 
the  locality  of  physical  properties  of  extended  systems  and 
provides  general  support  for  various  local  methods  of  calcu- 
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la t ion. 
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1.  Introduction 

This  note  deals  with  a particle  described  by  the  single  particle 
Hamiltonian 

H = j p2  + V(r)  . (1.1) 

where  V(r)  is  defined  for  all  r inside  a large  region  ft  which  eventually 
becomes  infinite.  The  eigenfunctions  and  eigenvalues  of  H,  <J>V(r)  and  EV, 
are  defined  by  the  Schroedlnger  equation 

H *v(r)  - Ev  i|iv(r)  . (1.2) 

the  boundary  condition 

r)  ■ 0 on  boundary  of  ft,  (1.3) 

and  the  normalization  condition 

f |<(/V(r)|2  dr  - 1 . (1.4) 

Of  course  a change  of  V(r)  in  the  vicinity  of  a point  rQ  in  general  changes 
the  eigenfunctions  ij/v(r)  at  all  points  r.  This  change  can  be  considerable 
even  if  the  distance  between  r and  ro  becomes  large.  For  example,  if  we 
start  with  V(r)  = 0 and  ft  a large  box  (-L  < x,y,s  < +L) , and  subsequently 
Introduce  a one  dimensional  localised  barrier 


V(r)  2 


-d  < x < d 
otherwise 


(1.5) 


the  eigenfunctions  ijiV(r)  are  changed  by  comparable  amounts  for  all  values 
of  r.  This  can  be  regarded  as  a kind  of  "action-at-a-distance"  of  the 
potential  V(r)  on  the  wavefunctiona  ipv(r). 


Let  us  now  define  the  quantity 


where  u is  an  imaginary  "time"  and  Pr(u)  denotes  path  Integration  over 


i r 


all  pacha  acarclng  and  ending  it  r » 0. 

Let  ua  divide  ft  Into  an  Interior  region 


ft. 


Inc 


and  an  exterior  region 


ext 


lx| , | y | . | s | $ 


lx| , |yl , |sl 


Wa  say  chan  write  0(3)  aa  Che  sum  of  two  terms. 


0(8)  - 0^(6)  + 0^(8)  , 


(2.2) 


(2.3) 


(2.4) 


where  0^(6)  comprises  all  contrlbutlona  to  (2.1)  coming  from  patha  entirely 

lnalde  & , and  p (3)  repreaanta  the  contrlbutlona  from  pacha  extending 

lnt  ext 

outaida  ft  . Wa  may  nota  that  P^nt(8)  la  the  LPF  for  a ayatem  with  potential 

V(r)  lnalde  ft,  _ and  Infinite  walls  at  the  boundary  of  ft, 

lnt  • int 

Now  consider  a comparison  potential,  V'(r),  identical  to  V(r)  Inside 

ft,  but  different  In  ft  . The  LPF  for  this  system  can  be  written  aa 
lnt  ext 


P'(8)  " Pint(6)  + 0,ext(B)  * 


(2.5) 

where  the  interior  contribution  la,  by  the  construction  (2.1)  the  same  aa 
for  0(8),  Eq.  (2.4). 

Thus  the  difference  of  the  two  LPFs  is  given  by 

Ap(8)  = O' (B)  - P(6)  - 0*^(6)  - Put(8)  . (2.6) 

Now  let  us  assume  that  both  V(r)  and  V'(r)  are  bounded  below 

V(r)  , V' (r)  > Valn  . (2.7) 


Then  clearly,  from  (2.1)  , 


• 8 

°ext(S)*  p’axt(S)  < J*3*  H <T  P(u)  + Vam)du|Pext  E(u) 


e“^vmin 


po,ext^  * 


(2.8) 


where  Pext  r(u)  denotes  integration  over  all  paths  which  extend  into 


a«t  and 


,.8 


po.«xc(6) 5 H |*j  i r<“>  <2-’> 


corresponds  to  a free  particle,  V(r)  = 0.  Next  let  us  assume  that  V(r) 
(but  not  necessarily  V’(r))  is  bounded  above. 


V(r)  < V. 


Than  clearly 


P(6)  > e~6Vn*x  po(g) 


(2.10) 


(2.11) 


where  PQ(S)  is  the  local  partition  function  for  a free  particle. 


/■ 

'0  • 

i/2 


PQ(8)  = J exp  J -J  j r*(u)du|pr(u) 

- 


(2.12) 


Combining  Eqs.  (2.6),  (2.8),  (2.11)  and  (2.12)  gives 


l&  P(8)  I ^ 2 e^^max  ~ vmin)  1 

P(8)  P„ 


Po,.xC<S> 


w 


(2.13) 


Lastly  we  calculate  0 (8),  the  contribution  to  P (8)  of  paths 

O i ®Xl  O 

extending  beyond  ^inC*  We  write 


Pn  ' MB)  - pn  . „(8) 

o,ext  o o.int 


(2. 1A) 
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and  recall  that  0Q  lnt(8)  is  the  LPF  for  free  electrons  In  the  box  (2.2) 
with  Infinite  walls.  The  cosine  eigenfunctions  and  eigenvalues  are 


“ ~T7T  cos 


3T(Vt+  -|)  *| 

ff(v2+  y)  y 

cos 

ff(v3+  -|)  z 

' I COS 

a J 

a 

a 

,v  1 rr2  f.  1 

2 7 I (vi  + 2 


) + (v2  + j)  + (v3  + 


(2.15) 

(2.16) 


where  v 5(Vj,  v2,  v3).  Eigenfunctions  which  are  odd  In  x,  y or  z do  not 
contribute  to  the  LPF  at  r ■ 0.  Therefore  we  have 


Oo.ln, «>  = *"  2 

v1,v2,Vj-0 


f 1 2 1 2 1 21 

2l*  {(M  2>  + <vz+  y>  + (vj+  y>  j 


2.  3 


- a-J  j£e_0&(y+  ] 

Using  Poisson's  animation  formula^  we  can  rewrite  the  sum  over  u as 

a_1£  e-6?!2^  2?  - (2ir6)"li(l-2  Z (-l)n+1  e ~ ) 


(2.17) 


so  that 


2nV 


PQ  .^.(B)  - (2TTB)"3/2{l-(l  - 2 D (-l)n+1  e'  6 )3  } 

O l "Xt  • 

n«l 


(2.18) 


(2.19) 


For  large  a only  the  leading  term  In  the  sum  needs  to  be  retained  giving 

-3/2  -iil 

Po,ext(0)  - (2ff8)  6 • 8 • (2-20) 

It  can  be  verified  that  the  right  hand  side  of  (2.20)  constitutes  an  upper 


bound  on  p _(B)  for  all  a, 
o, ext  1 ■ 


6 


°0,e*t(3)  4 <21t6)"T  6 6 


(2.21) 


Substitution  of  (2.12)  and  (2.21)  Into  (2.13)  gives  our  final  result 


1AP(8)I  ,,  e8(Vmax-Vmin)  e'  3 

P ( 6)  * e 


(2.22) 


In  the  Appendix,  a similar  result  is  proved  for  the  density  matrix  p(r ,r' ; 8) . 


Discussion 

The  key  feature  of  Eq.  (2.22)  is  the  Gaussian  decrease  of  Ap(0, 8)/p(0, 8) 

with  the  size  2a  of  the  cube  which  separates  the  interior  and  exterior 

regions.  Thus  we  see  that  any  change  of  the  potential  V(r)  outside  of  ft.  , 

- tnt 

subject  only  to  the  lower  bound  limitation  (2.7),  will  have  a negligible 
effect  on  p(0,S),  provided  that  a is  large  enough.  In  particular,  the  potential 
outside  of  ftint  may  be  replaced  by  zero  or  ftj may  be  enclosed  by  an  infinite 
barrier,  without  having  an  appreciable  effect  on  p(0,8). 

Since  P(r,8)  is  the  Laplace  transform  of  the  local  density  of  states 
n(r,E)  (See  Eqs.  (1.8)  and  (1.10))  the  "locality"  of  P(r,8)  implies  in  a 
general  way  the  locality  of  n(r,E).  Thus  the  theorem  (2.22)  which  we  have 
derived  provides  a general  support  for  the  physical  fact  that,  in  an  extended 
system  physical  properties  near  a point  r depend  significantly  only  on  the 
potential  near  r . It  also  lends  support  for  various  local  methods  of  calcu- 
lation, such  as  cluster  methods^3  moment  methods^*  ?\  etc.,  which  have 
been  extensively  used  in  recent  years  to  determine  local  properties  of  ex- 
tended systems. 

However,  the  Gaussian  decay  of  &p(r,8)  as  function  of  a,  by  no  means 
implies  that  the  inverse  Laplace  transform  of  P(r,8),  the  local  density  of 


I 


r 
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APPEND  LX 

In  Che  following,  we  prove  a localization  theorem  for  the  Laplace  transform 
of  the  density  matrix, 

p(r,r';t3)  = / n(r.r';E)  e_6EdE  , (A.l) 

where  n(r,r';E)  is  the  density  matrix  of  the  system  under  consideration, 

[ 

n(r,r';E)  -£>V*(r)  ^V’)  5(E-EV)  . (A.2) 

v 

The  Laplace  transform  p(r,r';S)  can  be  written  as  a path  incegral  similar  to 
Eq.  (2.1),  where  the  relevant  paths  are  those  starting  at  r and  ending  at  r'. 
Thus,  for  r.r'  inside  £2  p(r,r';3)  obeys  inequalities  analogous  to  Eqs.  (2.8) 

and  (2.13).  Consequently,  in  order  to  show  that  p(r,r';S)  is  localized,  it 
suffices  to  show  that  Po(r,r';S)  has  this  property,  where  p^  is  the  corresponding 
free-electron  quantity. 

The  contribution  of  paths,  entirelv  inside  the  cube  0,  ...  to  p can  be 

int  o 

evaluated  explicitly  by  using  (A.l)  and  (A.2)  and  the  appropriate  wave  functions. 
In  this  way  we  obtain 

0O  lnC  ■ J(x.x')  J(y.y')  J(z,z')  , (A. 3) 

where 

J(x.x')  -£tj(*(x)  ^Q(x’)  e_l3En  , (A. 4) 

tl 

and  Che  set  is  given  by 

_u 

(x)  - a cos (k  x)  ; k a - ir(n+  *i) 
in  in  2n 

*2n+i(x>  ’ ^ 3intk:n+ia)  : k*n+ia  " *(n+i) 

E - «s  k 2 ; n - 0,  1,  2 

n n 


- 


Separating  J into  even  and  odd  contributions,  and  using  Poisson's  summation 
formula,  one  can  show  that 


e 20 

J - r 1 + 

(2tt8)^ 


2E  e- 


8 cosh 


4an(x-x' 


_ 2a2 (2n-l) 2 f 

2 e 8 £ e 8 cosh  ■2-a(2n~1)  (x4* 

n-1  L 0 


(A. 5) 


For  a large  cube  of  edge  2a,  and  for  x and  x'  close  enough  to  the  origin  (in  such 
a way  that  |x|  + |x'|«  a),  the  series  in  (A. 5)  converge  very  rapidly  and  can  be 
replaced  by  the  corresponding  first  terms.  In  this  way  we  obtain 


J(i.x')  - J^x.x' ) jl  + 2«-  *f-  cosh[iiip^]-2e-  co.h[2‘^'>] 


(A. 6) 


where  is  the  corresponding  contribution  from  an  infinite  system,  i.e.. 


J_(x,x')  - (2ir8) 


, (x-x' ) 2 

-1*  e-^T^ 


(A.7) 


We  apply  now  Eqs.  (2.13),  (A. 3)  and  (A. 8)  to  obtain  the  following  localization 


theorem: 


i .‘1 
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